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Abstract 

In this second paper, we prove a necessity Theorem about the topological origin of phase tran- 
sitions. We consider physical systems described by smooth microscopic interaction potentials 
V]\r(q), among N degrees of freedom, and the associated family of configuration space sub manifolds 
{M v } V £hl, with My = {q € M. N \VN(q) < v}. On the basis of an analytic relationship between a 
suitably weighed sum of the Morse indexes of the manifolds {M v } ve ^. and thermodynamic entropy, 
the Theorem states that any possible unbound growth with N of one of the following derivatives of 
the configurational entropy S^\v) = (l/iV)log J Mv that is of \d k S ( '\v)/dv k \, for k = 3,4, 
can be entailed only by the weighed sum of Morse indexes. Since the unbound growth with N of 
one of these derivatives corresponds to the occurrence of a first or of a second order phase tran- 
sition, and since the variation of the Morse indexes of a manifold is in one-to-one correspondence 
with a change of its topology, the Main Theorem of the present paper states that a phase transition 
necessarily stems from a topological transition in configuration space. The proof of the Theorem 
given in the present paper cannot be done without Main Theorem of paper I. 

PACS numbers: 05.70.Fh; 05.20.-y; 02.40.-k 

Keywords: Statistical Mechanics, Phase Transitions, Topology 
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1. INTRODUCTION 



In Statistical Mechanics, a central task of the mathematical theory of phase 
transitions has been to prove the loss of differentiability of the pressure func- 
tion - or of other thermodynamic functions - with respect to temperature, 
or volume, or an external field. The first rigorous results of this kind are the 
exact solution of 2d Ising model due to Onsager 1|, and the Yang-Lee theorem 

showing that, despite the smoothness of the canonical and grand canonical 
partition functions respectively, in the iV — > oo limit also piecewise differen- 
tiability of pressure or other thermodynamic functions becomes possible. 

Another approach to the problem has considerably grown after the intro- 
duction of the concept of a Gibbs measure for infinite systems by Dobrushin, 
Lanford and Ruelle. In this framework, the phenomenon of phase transition 
is seen as the consequence of non-uniqueness of a Gibbs measure for a given 



type of interaction among the particles of a system 



Recently, it has been conjectured j], 



3,Q. 



, i| that the origin of the phase 
transitions singularities could be attributed to suitable topology changes within 
the family of equipotential hypersurfaces {Y> v = V^ 1 (f)} t , e K of configuration 
space. These level sets of Vn naturally foliate the support of the statistical 
measures (canonical or microcanonical) so that the mentioned topology change 
would induce a change of the measure itself at the transition point. In a few 
particular cases, the truth of this topological hypothesis has been given strong 

evidence: i) through the numerical computation of the Euler characteristic for 

ft 

the {S^j^gK of a two-dimensional lattice ip 4 model jjjj; ii) through the exact 
analytic computation of the Euler characteristic of {M v = V^ 1 ((— cxd, f])} ugR 
submanifolds of configuration space for two different models, the mean-field 
XY modelQ] and the fc-trigonometric model 

In the present paper we prove a necessity Theorem which implies that for a 
wide class of potentials (good Morse functions), a first or a second order phase 
transition can only be the consequence of a topology change of the submani- 
folds M v of configuration space, and this appears to be the truly primitive and 
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deep mathematical origin of the phase transition phenomena, at least for the 
mentioned class of potentials. 

The Theorem is enunciated as follows: 

Theorem Let Vn{qi, ■ ■ ■ , Qn) '■ R — > R, be a smooth, non-singular, finite- 
range potential. Denote by M v := V Iv 1 ((—oo,v}), dgI, the generic subman- 
ifold of configuration space bounded by £„. Let {q® e M N } ie [i,Af{v)] be the set 
of critical points of the potential, that is s.t. VViv(?)| w = ; and M{v) be 
the number of critical points up to the potential energy value v. Let T(qi , eq) 
be pseudo-cylindrical neighborhoods of the critical points, and fii(M v ) be the 
Morse indexes of M v , then there exist real numbers A(N,i,E ), gi and real 
smooth functions B(N,i,v,£o) such that the following equation for the micro- 
canonical configurational entropy S^ N \v) holds 



S { N \v) = ^log 



1 

N 



■ N 

/ d N q + ^A{N,i,E G ) gi Hi{M v 



+ Yl B(N,i(n),v-v^,£ ) 



n=l 



(details and appropriate definitions are given in Section^), and an unbound 
growth with N of one of the derivatives \d k S^\v)/dv k \, for k = 3,4 ; and 
thus the occurrence of a first or of a second order phase transition respectively, 
can be entailed only by the topological term J2iLoA(N,i,£o) gi fii(M v - eo ) + 
J2nS )+1 B(N,i(n),v-v v c (v \£ ). 

The above given expression for the entropy stems from a decomposition of 
the volume of a generic submanifold M v into two parts: the volume of the 
disjoint union of suitably defined neighborhoods of the critical points of the 
potential Vjv(g), and the volume of its complement. The latter is represented 
by the first term in square parentheses, and, after Main Theorem of paper I, 
it cannot entail unbounded growth with iV of up to the fourth derivative of 
the entropy. Only the second and third terms in square parentheses could do 
it. These two terms, representing the volume of the neighborhoods of critical 



points, are of topological meaning. 

Thus the proof of the present Theorem crucially relies on Main Theorem of 
paper I which is a first step toward proving the necessity of topology changes 
of configuration space submanifolds (either the level sets of the potential 
function, or the manifolds M v bounded by them) for the occurrence of phase 



transitions 
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kJ. Main Theorem of paper I is based on the assumption of 
the existence - at any number N of degrees of freedom - of an energy density 
interval [^0,%] free of critical values. Under this assumption, in paper I, we 
proved the uniform convergence in the thermodynamic limit of configurational 
entropy in the class of three times differentiable functions. However, since in 
general it is a very hard task to locate all the critical points of a given potential 
function, it is also very hard to ascertain whether Main Theorem of paper I 
applies to it or not. To overcome this limitation, in the present paper we prove 
a new Theorem, based on a less restrictive assumption (allowing the existence 
of critical points), which we could not prove without the Main Theorem of 
paper I. 

As already remarked above, the necessity Theorem proved in the present 
paper applies to a very broad class of systems: those described by finite range 
potentials which are good Morse functions. In fact, checking whether a given 
potential is a good Morse function or not is not difficult and amounts to con- 
trol whether the potential is smooth, bounded below and whether its Hessian 
is non-degenerate, that is free of vanishing eigenvalues; degeneracy typically 
occurs in presence of continuous symmetries and can be removed by arbi- 
trarily small and standard perturbations [14] which do not alter neither the 
microscopic dynamics produced by the gradients of the potential nor the ther- 
modynamics. 

2. BASIC DEFINITIONS 

For a physical system S of n particles confined in a bounded subset A d of 



l d , d = 1,2,3, and interacting through a real valued potential function Vj 



N 
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defined on (A d ) xn , with N = nd, the configurational microcanonical volume 
Q(v,N) is defined for any value v of the potential Vn as 

Q(v,N)=f dq 1 ...dq N 5[V N (q 1 ,...,q N ) -v]= [ ° , (1) 

</(A d )x™ Js„ ll vl/ Anl 

where da is a surface element of E„ := V^ 1 (f); in what follows Q(v,N) 

is also called structure integral. The norm ||VV/v|| is defined as ||VVjv|| = 

[Eii(^^v) 2 ] 1/2 - 

Now we can define the configurational thermodynamic functions to be used 
in this paper. 

Henceforth, according to the need for explicit reference to the iV-dependence 
of V, we shall use both V and Vn to denote the potential. 

Definition 1. Using the notation v = v/N for the value of the potential energy 
per particle, we introduce the following functions: 

- Configurational microcanonical entropy, relative to E^. For any JVgN 
and i)6K, 

S N (v) = S N (v; V N ) = 1 logQ(Nv, N) . 

- Configurational microcanonical entropy, relative to the volume bounded 
by E„. For any N e N and t)6l, 

S { N \v) ee S { N \v;V N ) = 1 log M(Nv,N) 

where 

da 



M(v,N)= [ dq 1 ...dq N e\y N (q 1 ,...,q N )-v]= f drj [ 

J(A d ) x " JO Vs r 

(2) 

with Q[-] the Heaviside step function; M(v,N) is the codimension-0 subset of 
configuration space enclosed by the equipotential hypersurface E„. The repre- 
sentation of M(v, N) given in the r.h.s. stems from the already mentioned 
co- area formula in r L 15]. 

Definition 2 (First and second order phase transitions). We say that a 
physical system S undergoes a phase transition if there exists a thermodynamic 
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function which - in the thermodynamic limit (N — > oo and vol(A d ) / N = const) 
- is only piecewise analytic. In particular, if the second-order derivative of 
the entropy S^\v) is discontinuous at some point v c , then we say that a 
first-order phase transition occurs. If the third-order derivative of the entropy 
sL\v) is discontinuous at some point v c , then we say that a second-order 
phase transition occurs. These definitions stem from the standard definitions 
of first and second order phase transitions as due to a discontinuity of the first 
or second derivatives of the Helmoltz free energy, respectively, and from the 
existing relationship - through a Legendre transform - between the Helmoltz 
free energy and the entropy (see Definition 1 in paper I). 

Definition 3 (Standard potential, fluid case). We say that an N degrees 
of freedom potential Vn is a standard potential for a fluid if it is of the form 

V N : B N C R N -> R 

n n 

V N (q) = £ *(||ft-£||) + 5>A(3) (3) 

i=£j=l i=l 

where Bn is a compact subset ofM N , N = nd, ^ is a real valued function of 
one variable such that additivity holds, and where Ua is any smoothed potential 
barrier to confine the particles in a finite volume A, that is 

ifqeA' 
U\(q) = < +oo if q E A c , complement in M. N 
I C°° function for g G A \ A' 

where A'c A and A 1 arbitrarily close to A C 1^, closed and bounded. Ua is a 
confining potential in a limited spatial volume with the additional property that 
given two limited d-dimensional regions of space, A 1 and A 2 , having in common 
ad— 1-dimensional boundary, Ua 1 + Ua 2 = ^Aiua 2 - By additivity we mean 
what follows. Consider two systems S\ and S 2 , having Ni = n±d and N 2 = n 2 d 
degrees of freedom, occuping volumes Af and A d , having potential energies v 1 
and v 2 , for any (q 1 , ...,q Nl ) e (Af) xni such that V Nl (qi, . . . , q Nl ) = v x , for 
any (q Nl +i, ...,qN 1 +N 2 ) G (A^) x ™ 2 such that V N , 2 (q Nl+1 , . . . , ?at 1+ at 2 ) = v 2 , for 
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(qi, . . . , q^+Ni) £ (Af) xni x (A d ) xri2 let V N (q 1 , . . . , qN 1 +N 2 ) = v be the potential 
energy v of the compound system S = S\ + S 2 which occupies the volume 
A d = Af U A2 and contains N = Ni + N 2 degrees of freedom. If 

vim + N 2 , A d U A|) = v^N,, A?) + v 2 (N 2 , A d 2 ) + v'(N u N 2 , A d , A d 2 ) (4) 

where v' stands for the interaction energy between Si and S 2 , and ifv'/vi — > 
and v' jv 2 — > /or AT — > 00 i/ien V)v additive. Moreover, at short distances \& 
mttst fre a repulsive potential so as to prevent the concentration of an arbitrary 
number of particles within small, finite volumes of any given size. 

Definition 4 (Standard potential, lattice case). We say that an N degrees 
of freedom potential Vn is a standard potential for a lattice if it is of the form 

V N : B N c R N -> R 

v N { q )= c anWi-Qj\\)+ E *(® (5) 

ijeXcN d iei cN d 

where is a compact subset ofR N . Denoting by ai,...,aa the lattice spac- 
ings, if i G N d , then (iiai, ... ,140,4) G A d . We denote by m the number of 
lattice sites in each spatial direction, by n = m d the total number of lattice 
sites, by D the number of degrees of freedom on each site. Thus q\ G R D for 
any i. The total number of degrees of freedom is N = m d D. Having two sys- 
tems made of N = m d D degrees of freedom, whose site indexes i^ and run 
over 1 < i\\ . . . , < m, and l<if\..., i^ < m, after gluing together the 
two systems through a common d — 1 dimensional boundary the new system 
has indexes i running over, for example, 1 < i± < 2m and 1 < i 2 , . . . , 14 < m. 
U 

v(N + N, A d U A*) = Vl (N, A d ) + v 2 (N, A d ) + v'(N, N, A d , A d 2 ) (6) 

where v' stands for the interaction energy between the two systems and if 
v 1 / v\ — > and v'/ v 2 — > for N — > 00 then Vn is additive. 

Definition 5 (Short-range potential). In defining a short-range potential, 
a distinction has to be made between lattice systems and fluid systems. Given 
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a standard potential Vn on a lattice, we say that it is a short-range potential if 
the coefficients are such that for any i,j G X C N d , Cjj = iff \i — j\ > c, 
with c is definitively constant for N — > oo. 

Given a standard potential Vn for a fluid system, we say that it is a short- 
range potential if there exist Rq > and e > such that for ||q|| > Rq it is 
I^GIqIDI < IIqII - ^ 6 '') where d = 1,2,3 is the spatial dimension. 

Definition 6 (Stable potential). We say that a potential Vn is stable 
if there exists B > such that 



for any N > and (qx,... 7 q N ) G (A d ) xn , or for q} G R D , i G X C N d , 
N = m d D, for lattices. 

Definition 7 (Confining potential). With the above definitions of standard 
potentials Vn, in the fluid case the potential is said to be confining in the sense 
that it contains which constrains the particles in a finite spatial volume, 
and in the lattice case the potential Vn contains an on-site potential such that 
- at finite energy - ||^|| is constrained in compact set of values. 

Remark 1 (Compactness of equipotential hypersurfaces). From the 
previous definition it follows that, for a confining potential, the equipotential 
hypersurfaces E„ are compact (because they are closed by definition and bounded 
in view of particle confinement). 

In view of formulating and proving the Main Theorem of the present paper, we 
have to define some neighborhoods, that we call "pseudo-cylindrical" , of critical 
points of a potential function Vn- Before defining these pseudo-cylindrical 
neighborhoods of critical points, let us remember the following basic result in 
Morse theory. 

Theorem. Let f be a smooth real valued function on a compact finite dimen- 
sional manifold M . Let a < b and suppose that the set 




(7) 



f-\[a, b])=M = {xe M\a < f(x) < b) 



(8) 
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is compact and contains no critical points of f , that is ||V/|| > C > with C 
a constant. Let y G (a, 6) . Then there exists a diffeomorphism 

a:(a,b)x f-\y)^ f-^b)] by (v, x) >- a(v, x). (9) 

Corollary. The manifolds f~ 1 {y), a < y <b, are all diffeomorphic. 

This result is based on the existence of a one-parameter group of diffeomor- 
phism 

a v : M ^ M by x >-» <r(u, re) (10) 

associated with the vector field X = V/(x)/||V/(x)|| 2 with v >— » a(t',x) a 
solution of the differential equation on M 

da(v,x) Vf[a(v,x)) 

dv \\Vf[*(v,x)W> a[ ^ X) X - 

a(v,x) is defined for all v G R and x G M. Details can be found in standard 
references as 
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Applied to the configuration space M, if the function / is identified with 
the potential Vn, then in the absence of critical points of V in the interval 
(vq,vi) the hypersurfaces £„ = V / ^ 1 (f), v G (vq,V\), are all diffeomorphic. 

Definition 8 (Pseudo-cylindrical neighborhoods). Let E Vc fre a critical 
level set of Vn, that is a level set containing at least one critical point of Vn- 
Around any critical point qi\ consider the set of points ^{qi , p; v c ) C S„ c at 
a distance equal to p > from q£ , that is q G j(qc , p; fc) =^ ^(<Z — 9c ) = p ; 
where d(-,-) is the distance measured through the metric induced on S„ c fry i/ie 
euclidean metric of the immersion space, and p is such that p < | mhijj d(qf' ' — 
g c ), J ^afre/ i/ie critical points on the given critical level set. Moreover, set 
the thickness of all the pseudo-cylinders equal to Eq = minj e fq(v J c +1 — v{). After 
Sard Theorem, both p and Eq are finite because, at finite dimension, there is a 
finite number of isolated critical points and, consequently, a finite number of 
critical values. We define a pseudo-cylindrical neighborhood T(qc ,£q) C M of 
qc as the open subset of M bounded by the following set of points. By mapping 
j(qc\p;v c ) from S t , c to £(„ c+£o ), and from T, Vc to Y> {Vc _ £o) , through the flow 
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generated by the vector field X = VVjv(g)/|| VVAr(g)|| 2 , we obtain the walls of 
T(qc\e ), which are transverse to the T, v , and then we close the neighborhood 
with the pieces o/E(„ c+eo ) and E(„ c _ eo ) bounded by the images 7((?c , p; (v c + £o)) 
and 7(?c , P', (v c — eo)) of ^(qc\ p; v c ) through a(v, x), respectively. 

Lemma 1 (Generalization of Corollary 1, paper I). Let Vn be a standard, 
smooth, confining, short-range potential bounded from below (Definitions^^ 

V N : B N c R N -> R 

with Vjv given by Definition^ (fluid case), or by Definition^ (lattice case). 

Let {E„}„ g R be the family ofN — 1 dimensional hypersurf aces := V^ l {v), 
v G R, of W N . Let {M„}„ e t 6e i/ie family of N dimensional subsets M v := 
oo, v\), v G R, ofR N . Let {M v } ve ^ be the family of N dimensional 
subsets M v := M„ \ r(g£\ e), v G R, 0/ R^, w/*ere r(gi <} ,e) are tfie 

pseudo-cylindrical neighborhoods of the critical points qi^ o/V/v(o) contained in 
M v and J\f(v) is the number of critical points in M v . Let {E„}„gR be the family 
ofN — 1 dimensional subsets ofM N defined as E„ := E. u \ U^f^[r(gc , e) H E v ] . 
Let Wo = Vo/N,vi = Vi/N G R, ^0 < ^1 onc ^ ^ — v c /N be the only critical 
value ofVN in the interval 1 \ = [vq,vi], and let T*(qc ,£*), with q£ G V^ 1 (f c ) 
and e* such that e* > max(t>i — v c , v c — Vq). The following two statements hold: 

a) for any v, v' G [vo, vi] it is 

E m is C°° — diffeomorphic to E m r ; 

b) putting M(v,N) = vol(M v ), the quantities [d~M(v, N)/dv]/M(v, N) and 
(d k /dv k ){[dM(v,N)/dv]/M(v,N)}, k = 1,2,3, are uniformly bounded in N 
in the interval [vq,vi]. 

Proof. For what concerns point (a), we note that the flow associated with 
the C°° vector field X = VVjv(g)/|| VVAr(g)|| 2 is well defined at any point 
q G M Vl \ M vo . Thus the set M V1 \ M vo is diffeomorphic to the non-critical 
neck dM V0 x [t> ,t>i]. Then, after the "non-critical neck theorem" ^?|], for 
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any v, v' G h = [vq,Vi] it is « Incidentally, this entails also 

Mnv ~ Mnv' for any v,v' G [t?o,£>i]. 

Now let us consider point (b). Define Sn(v) = j^log[fl(Nv, N)}, where 
Q(Nv, N) = vol(E B jy)]> having proved the statement (a), we can apply Lemma 
4 of paper I which entails that 



sup ISVtv)! < oo and sup 



d k S 



dv 



V] 



< oo, k = 1,2,3,4. 



Whence, after Lemma 3 of paper I, it follows Soo(v) = linijv^oo Sn(v) G C 3 (/^). 
The next step is to prove that also sU(v) G C 3 (!$)■, where 



S ( J(v) := lim S ( -\v) = lim 1 log[vol(M^)], 

because, after Lemmas 3 and 4 of paper I, this entails the truth of statement 
(b). Let us begin by considering the microcanonical configurational inverse 
temperature. From its definition (5^{v) = dS-fr /dv one obtains Pn(v) = 
Q(Nv, N)/M(Nv, N). The function f3^{v) is well known to be intensive and 
well defined also in the thermodynamic limit, at least for extensive potential 
energy functions. Then we work out a representation of /3n{v) in the form 
of a microcanonical average of a suitable function. To this purpose we derive 

nn 

fl(Nv, N) with respect to v by means of Federer's derivation formula |15l l20j|. 
and then we integrate it. Federer's derivation formula (see Lemma 5 of paper 
I) states that 

>-i l,v "WA • (12) 

where A stands for the operator 

W 



A(.) = V 

Then we can write 



W J W 



n(Nv,N) = f drj f ||W|| A ( 

JO J'Er, V 



da 



M v 



WW J IIWI 

wi 



dji \\W\\ A I -fj^r— ) (13) 
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where dfi = d N q, so that we finally obtain 



Pn(v) 



-i -l 



M v 



dfx ||W|| A — — = ( || W|| A 



w\W\\j V 11 \\\vv\ u/Mn 



(14) 

which holds for v G [0,Uo). An important remark is in order. We have used 
Federer's derivation formula apparently ignoring that it applies in the absence 
of critical points of the potential function. However, if the potential V is a 
good Morse function (not a very restrictive condition at all) we know, after 
Sard theorem 17[, that the ensemble of critical values, here of the potential, 
is a point set. Therefore, any finite interval of values of the potential is the 
union of a finite number of open intervals where no critical value is present, and 
correspondingly no critical point on the exists. On all these open sets, free 
of critical points, Federer's derivation formula can be legally applied. Moreover, 
the results found by applying Federer's formula on each open interval free of 
critical values of V can be regularly glued together because of the existence of 
the thermodynamic limit of /3n(v). 

Let us now consider Cl(Nv,N) for v G [uo,ui]. As all the hypersurfaces 
S c labeled by v G [^o; Ui] are diffeomorphic, we can use Federer's derivation 
formula to obtain an expression for Q(Nv, N) similar to that given in Eq. (fH 
for Q(Nv,N), that is 

TI(Nv,N) = jf d V J_ ||W|| ^(p^Ji) +V(mo,N) 

dfi ||W|| A ( j-^—r J + / dfx A 



M V \M VQ \\\VV\\J J Mvq ' \I|VV| 

dp \\VV\\ A (t^) (15) 
'Af„\r* Vll vv/ ll/ 

where T* stands for the union of all the pseudo-cylindrical neighborhoods of 

the critical points of V in the interval [vq, v±\. Then we consider the restriction 

(3 N {v) of the function Pn(v) to the subset M v \ T*; from 

NK 1 M{Nv,N) v ; 



13 



we get 



M V \T* 

W\\ A 



m v \t* 



W 



dp \\W\\ A 



M„\r* 



W 



(17) 



By comparing Eq. fHj) with Eq. (fTY|) . we see that also P^(y) has to be intensive 
up to the iV — > oo limit, like 0n(v). In fact the excision of the set T* out of M v , 
no matter how the measure of T* depends on N, cannot change the intensive 
character of (3 N {v). The relationship among Sn(v), S n (v) and (3 N (v) is given 
by the logarithm of both sides of (fTB|) 

1 logO^iV, N) = ^ \ogM(vN, N) + 1 log^ N (v) . 

whence, using lini7v->oo jjlog(3 N (v) = 0, we obtain (v) = Soo(v) and thus 

st\v)eC\h). 

-(-) 



Finally, (v) e C (la) entails 



sup 



S ( N \v) 



< oo and sup 



< oo, k = 1,2,3,4. 



so that, resorting to Lemma 3 of paper I, the truth of statement (b) follows. □ 



3. MAIN THEOREM 



In this Section we prove the following 

Theorem 1 (Entropy and Topology). Let Vjv(<7i, • • • , Qn) '■ — >• M, be a 

smooth, non-singular, finite-range potential. Denote by M v := V^ 1 ((— oo, v ]), 
!)6l, the generic submanifold of configuration space bounded by £„. 

Let {q® G R N } ie [i,Af(v)] be the set of critical points of the potential, that is 
s.t. VVat(9)| w = 0, and J\f(v) be the number of critical points up to the 
potential energy value v. Denote by v = v/N the potential energy density, 
and assume that for any given interval [vq,vi] the number of critical values 
v{ contained in it is at most a linearly growing function of N . Let r(gc , £o) 
be the pseudo- cylindrical neighborhood of the critical point qi l \ and /ii(M v ) be 
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the Morse indexes of M v , then there exist real numbers A(N,i,Eo), g% and real 
smooth functions v >— > B(N,i,v,e ) such that the following equation for the 
microcanonical configurational entropy S^\v) holds 

N 

MAT 4 c „\ r,. 11 A A/, 

l V—£QJ 



S ( -\v) = ^lof 



N 



d N q + J2 A (N,i,e ) ft fii(M v 



i=0 



+ E 5 ^ 



i(n),v-v v J v \e ) 



n=l 



18) 



where v{y) = max{j\v J c < v}, and an unbound growth with N of one of the 
derivatives \d k S^\v)/dv k \, for k = 3,4, and thus the occurrence of a first 
or of a second order phase transition respectively, can be entailed only by the 



topological term Y^f =0 A(N,i, e Q ) g { /ij(M„_ £0 ) + J2n=i B(N,i(n),v-Vc (v) ,e ). 



The proof of formula (|T%j) is worked out constructively. This formula re- 
lates thermodynamic entropy, defined in the microcanonical configurational 
ensemble, with quantities of topological meaning (the Morse indexes) of the 
configuration space submanifolds M v = Vy 1 ((— oo, v]) = {q = (qi, . . . , qjy) G 
R N \V N (q) < v}. 

After Morse theory, topology changes of the manifolds M v can be put in 
one-to-one correspondence with the existence of critical points of the potential 
function Vjsf(qi, . . . , q^)- A point q c is a critical point if VVAr(g)|q =gc = 0. The 
potential energy value v c = V^{q c ) is said to be a critical value for the potential 
function. Passing a critical value v c , the manifolds M v change topology. Within 
the framework of Morse theory, if the potential Vjv is a good Morse function, 
that is a regular function bounded below and with non-degenerate Hessian 
(that is the Hessian has no vanishing eigenvalue), then topology changes occur 
through the attachment of handles in the neighborhoods of the critical points. 
Therefore, in order to establish the relationship between entropy and configu- 
ration space topology, we have to unfold the contribution given to the volume 
of M v by suitably defined neighborhoods of all the critical points contained 
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in M v because it is within these neighborhoods that the relevant information 
about topology is contained. 

This result is made possible by the idea of exploiting the existence of the 
so-called Morse chart in the neighborhood of any nondegenerate critical point 
of the potential function Vn- In fact, the Morse chart allows to represent 
the local analytic form of the equipotential hypersurfaces in an universal form 
independent of the potential energy value at the critical point, and only de- 
pendent upon the index of the critical point (equal to the number of negative 
eigenvalues of the Hessian of the potential) and, obviously, upon the dimension 
iV of configuration space. Hence the possibility of a formal computation of the 
contribution of the neighborhoods of all the critical points to the volume of 
M v as a function of v. 



Proof. Let us consider the definition of the configurational microcanonical 
entropy S N \v), already given in Eq.(J2J), 

S { N \v) = ^logM(v,N) , (19) 

with 



V 



M{v,N)= I d N q = dr) d N q S[V N {q) - rj\ 

<V N (q)<v JO J(A d ) xn 

„ *• L • (20) 

where we have set equal to zero the minimum of the potential. Let 
be the family of all the critical level sets (in general not manifolds) of the 
potential, that is the constant potential energy hypersurfaces that contain at 
least one critical point q¥\ where VVjv(5)| (»} = 0. For a potential which is 
a good Morse function, after the Sard Theorem (see Corollary 2 at p. 200 of 
Ref.Q), at any finite dimension N, and below any finite upper bound of the 
potential energy, the number of critical points in configuration space and thus 
also the set of critical values {f^} je N, are finite, isolated and such that v{ < v* 
if j < k, so that any energy interval Vq < v < V\ is the union of a finite number 
of open intervals free of critical points. 
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In order to split the integration on M v into two parts: the integration on 
the union of the neighborhoods of all the critical points contained in M v and 
the integration on its complement in M v , we have defined for each critical 
point qi its pseudo-cylindrical neighborhood T(qc , £o)j £ o is the thickness - 
in potential energy - of the neighborhood. The assumption that the number 
of critical values v{ is at most linearly growing with iV entails, together with 
Sard theorem, that Eq is finite. 

Let us now split the integration on M v into the integration on M v f][J i T(qc\e ) 
and on its complement M v \ [J i T(qc\eo). We have 

d N q = [ d N q + [ d N q , (21) 

M v Jm v \(J?J? + ^ r(# ,e ) J M v n[jt { i +£0) r (^> .eo) 

where N(v ) is the number of critical points of Vn(q) up to the level v. We can 
equivalently write 

Af cl (v+e ) Nip 

vol(M v )= d N q+ Yl E/ , , ( 22 ) 

where Af c i(v) is the number of critical levels E^ such that u£ < u, and A/^, is the 
number of critical points on the critical hypersurface YP and where we have 
changed the notation of the pseudo-cylindrical neighborhoods to Tj(qi m \ e ) 
labelling with j the level set to which it belongs and numbering with m the 
critical points on the j — th level set. Notice that J\f(v) = Y^!j=x ^ -^cp- 

Then we use the co-area formula in the r.h.s. of Eq. (J2(J|) to rewrite Eq. (1221) ; 
a distinction is necessary between two cases for E„ = dM v : its label v is closer 
than Eq to a critical level or not; thus we obtain 

vol(M„) = / d N q 

J M v \[jt ( i +£0) ni ( c\e ) 
JV c i(_v+e ) Ml v V 3 c+£Q 

+ E E / ■ dr > / , , W?) " $ (23) 

J=1 m=1 Jv J c -e JTj(ql m ',e ) 

when t> > t>c^' ) + Eq and i> < v"^ +1 — Eq, where v(v) is such that < v < 
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; whereas 



vol(Af„) = / d N q 

iA/AU^ ( r eo) r( g «,e ) 

Af cl (v)-i Nip i +eo 

+ E E /. ^ / , , 

„-_i ™_i Jvk—en JTilqi™' ,en) 



j=l m=l Jv t- £ o J~Tj(<li Vo) 

+ E /„ d W , , AWg)-il- (24) 

m =l 17 «c -eo • / r„ (t ,)(g^ ' ,e ) 

when t> c — sq < v or u < v c + sq. 

Near to any critical point, a second order power series expansion of V(q) 
reads 

vS\q) = V N (q c ) + \ £ ^ W - ql) W - 4) 

For sufficiently small e , the integrals J r ^ £ ^ d N q 5[Vjy(g) — r/] can be replaced 
with arbitrary precision by j r ^ ^ ePg 5[Vjy 2 (?) —rj\. Moreover, if Vjv(g) is a 
good Morse function, then a coordinate transformation exists to the so-called 

k N 

V^(x) = V N (q c )-J2^i+ E 
i=l z=fc+i 

where /c is the Morse index of q c . Using Morse chart we have 



<Tg -v}= d N x | det J| 6[V^(x) - rj] (25) 

rj(g c ,eo) Jrj(q c ,eo) 

where J is the Jacobian of the coordinate transformation. 

Using Morse coordinates inside the pseudo-cylinder Tj(qi m \ e ) around the 
critical point gc m \ we see that each part of an hypersurface T, v fl Tj(qi m \e ) 
is a quadric 

fc m iV 

e=»7-^=-x;^+ E * 2 =-m 2 +m 2 , (26) 

Z=l l=k m +l 

where the Morse index of qi^ is denoted by k m , so that | X \ 2 = Y^i=i x m anc ^ 
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I Y \ 2 = YliLk +1 x f- Thus we rewrite the r.h.s. of Eq. (J25j) as 

|detJ| / d N x 5(-|X| 2 + \Y\ 2 - f ) 

= | det J| 

I dQ km ~ 1 dQ N - km - 1 d\X\d\Y\\X\ km - 1 \Y\ N - km - 1 5(-\X\ 2 + \Y\ 2 - £) 

(27) 

where dfl r is the solid angle element in r dimension, whose integration yields 
the surface C r of the r-dimensional hypersphere of unit radius. Putting z = 
\X\ 2 and integrating on the angular coordinates we get 

1 r<*(£,r) rP(€,r) 

-\ det J\C N _ km ^C km ^ dz / d\Y\\Y\ N - k -- 1 z^- 2 ^ 2 5{-z+\Y\ 2 -^) 

1 JO Jy^ 

(28) 



where a(£,r) = ( v^ 2 + 4r 2 - 0/2 and /3(f,r) = ^ (^/C 2 + 4r 2 + 0/2. These 
expressions stem from the definition of Tj(q c , Eq) whose boundaries have to be 
orthogonal to the potential level sets described by Eq. fpjj) . These boundaries 
are given by the equation |-X"||y| = r. Putting y = \Y\, from Eq. (|2"Hj) when 
£ > we obtain 

i J jm CW^-i / dy y N - k ^ (y 2 - O^' 2 (29) 
and when ( < we obtain 

-J ]m C^-iO^-i / dy y N ~ k ^ (y 2 - (fcm ~ 2)/2 (30) 



where CV-fc m -i an d C^ m are surfaces of hyper spheres of unit radii, that is 
C n = 27r"/ 2 /(n/2 - 1)! (for n even) and C n = 2("+ 1 VV n - 1 )/ 2 /(n - 2)!! (for 
n odd); Jj m stands for the numerical absolute value of the determinant of J 
computed at the critical level v{ and at the critical point qc ■ By defining 

F+(t,k m ,N) = / dy y"-^" 1 (y 2 - {km ~ m (31) 

and 

r/Jtt.r) 

F.(t,k m ,N) = / dy y^-- 1 (y 2 - {km ~ 2)/2 (32) 
Jo 
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we can now write 



vol(M„) = / d N q 



M c i(v+s ) Mi 



+ E Y<2 CN ^- lCkm - lJ i m dtF{t,k m ,N). (33) 

j=l 771=1 ^- £ 

when t> > Uc + 6 and t> < — £o> where = max{j|t^ < t>}, or 

vol(M v ) = / d N q 
N cl {v)-1 Nl v 1 

+ E Y,o CN - km - lCkm - lJjm diF{i,k m ,N) 

j=l m =l •'-eo 

+ E o^ C w J ]™ d£F(Z,k m ,N) (34) 

m=l J - £ o 

with j = J\f c i(v), when t>c^' )+1 — £o < f or t> < Uc^+£q. In Eqs. (j3^j) and (J54*j) we 
have put £° £o d£ ^(e, Av™, ^) = f_° £o ^ N) + ^ d£ F + {£, k m , N). 

Notice that it is J\f(v) = Y^f=o fH{M v ) , where Hi(M v ) are the multiplicities 
of the critical points of index i (there are at most N +1 values for the indexes 
of critical points at dimension N) below the energy value v. 

Therefore, we can rearrange the double summation in Eqs.()33|). (J54"jl by 
expressing it as a double summation on all the possible values of the Morse 
indexes and on the number of critical points for each value of the Morse index, 
that is 

(35) 

i=0 k=l 

where, since the integrals in Eqs.(|33|l. (J54*j) are independent of the index j, we 
have defined a set of positive coefficients A(N, i, Eq) as 



1 



A(N, i, e„) = -Cjv-f-iCi / di F{t, i, N) . (36) 

We remark that the term ()35|). being a function of the Morse indexes Hi(M v ), 
has topological meaning. In order to make clearer this topological meaning, 
we rewrite ()35|) in the equivalent form of a weighed sum of Morse indexes as 
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follows. From the set of positive numbers Jj(i t k) m (i,k) 



we define 



9i 



A M\ E ' / '"- 



Hi(M v 



k)m(i,k) 



k=l 



and rewrite the second term of the r.h.s. of Eq. (}3*3*j) as 



A' 



^A(N,i,e ) 9i m(M v ). 



(37) 



(3? 



Equation (J35)) is here recast in a form which makes somewhat clearer its 
dependence on Morse indexes. Moreover, we introduce also the coefficients 



B(N, i, v - <M e ) = ^C^QJ-^ 



1 



*(0) 



deF(e,z,iv), (39) 



so 



where fc(z) stems from k) = j, such that for v = Vc it is B(N, i, 0, Eq) = 0, 
and for v — v" = £o it is B(N, i, eq, eq) = A(N, i, £o)gi- 

For the purposes of the present proof, we are not concerned about the 
complication of the coefficients A(N, i, Eq) and B(N, i,v — v v c , Eq) because all 
what we need, in order to make the link between configuration space topology 
and thermodynamics, is that the second term in the volume splitting in Eq. (J21J) 
can be written in the form ()38|) . In fact, now we can write the entropy per 
degree of freedom as 



S { N \v) 



ilogM(« J JV) = llogj^ d N q 



(40) 



— log 



d N q+ / d N q 

a/AU^ +£0) r(#, £0 ) J M v n[jt ( i +£0) r(g£U) 



N 



log 



A? 



i=0 



i ^ u(v)+l u(v) 

when v < v c — Eq or v > v c + Eq, or 



4 _> («) = x l0i 



N 



M„\U^ +£o) r( g «, e o) 



rf Jv g + ^A(AT,i,e ) gifii(M v - 



so J 



i=0 



v (tj+£q) + 1 



+ E ^ 



.^n),u-< w ,eo) 



71=1 



(41) 
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I ^ u(v)+l , v(v) . 

when v > v c — £o or v < v c + sq. 

The equation above links thermodynamic entropy with the Morse indexes 
of the configuration space submanifolds M v , that is with their topology. In 
fact, according to Bott's "critical-neck theorem" (21 [, any change with v of any 
index /Zj(M„), i = 0, . . . , N, which can only be due to the crossing of a critical 
level, is associated with a topology change of the M v . 

Conversely, any topology change, in the sense of a loss of diffeomorphicity, 
occurring to the M v when v is varied, is signaled by one or more changes of 
the Morse indexes Hi(M v ) because, after the "non-critical neck theorem" [17] ]. 
this has to be the consequence of the crossing of a critical level. 

We remark that the two terms ^2 ^i9i^i an d Yli &i in Eq. (j41j) stem from the 
same term in the volume splitting (|2*Tj) (the union of neighborhoods of critical 
points), so that they will both participate in producing the development of 
singularities proper to a phase transition (for more details about the role of 
these two terms, see Remark EJ). 

Let us now show that the coefficients B(N, i,v — i>c , £0) are smooth func- 
tions of v at any finite N. 

Noting that d k B(N,i,v,e )/dv k = S k ~^F(^, k m , N)/d^ k ~ l \ we fo cus on 
the smoothness of F±(£,h m ,N). Let us consider F + (£,k m , N). There are 
two cases: k m is even so that (k m — 2)/2 is an integer; k m is odd so that 
{k m — 2)/2 = n + 1/2 with n e { — 1} UN. In the first case, by itera- 
tively applying the derivation formula ^ J^a) dxf(x,a) = f(ip(a), a) — 
d< ^^ f(4>(d), a) + dx d ^Q a a ^ to Eq. (jSU, one immediately realizes that 
F(£, k m , N) is smooth. 

In the second case, Eq.f^TJ) is rewritten as 

I(t,n,N) = / dy y"-'-- 1 (y 2 - n+ ~ 2 (42) 
which, after integration by parts, yields the recursion formula 



(y2 _ C)n+I I 



2 + 2n + N 



+ jt^Tn i( ^ n - 2) - (43) 
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For this recursion formula there are two possible initial conditions. The first 
initial condition is obtained by direct integration of Z(£, n, N) for N = k m + 2 
and n > 0, and is 



l^,n,N = k m + 2) 



(y 2 - r 



3 + 2n 



(44) 



the second initial condition for the recursion ()43|) is obtained by working out 
Z(£, n, AT) for N = k m + 1 and n > 0. From the above definition of Z(£, n, AT) 
we get 



z(e,n,Ar = fc m + i; 



(v 2 - n+l 



2(n + l)y 



(45) 



where the integral in the r.h.s. - that we denote by I(n) - can be solved by 
introducing the double recursion 



Kn) 



2(y 2 -ey 



(An + 3)r 



3£ 



4n + 3 



-#(n) 



2/=/3(£,r) 

H(n) = I(n-l)-£H(n-l) 



(46) 



which leads to 

H(n + 1) 

together with 



2(y 2 -O r 



(An + 3)y 3 



3£ 



An + 3 



+ Z)H(n) 



(47) 



H(l) 



(3? " 4) v ^~^ + log(y + v 7 ^ 3 ^) 



-log(v^). (48) 

From Eq. (}4*%)) and by Eq. (|4T|) the system of double recursion (J4T)|) is 
solved. From this we get the required second initial condition to be inserted in 
Eq. (|43p. Now, notice that all the terms that are either powers of (y 2 — £) or 



simply powers of y, are to be computed at y — /3(£, r) = y (a/£ 2 + 4r 2 + 0/2, 
so that they are never vanishing functions of £, as a consequence these terms are 
infinitely many times differentiable with respect to £. The apparent singularity 
of i?(l) for £ = is cured in Eq. (j4T)j) where it enters multiplied by £. We 
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remark that the recursion formula given in Eq. (|43|) holds also for n = —1,0 
which correspond to Morse indexes k m = 1, 2 respectively; the initial conditions 
for the recursion are to be explicitly computed from Eq. (|42|) by substituting 
n = — 1 or n = 0. Also for these special values of n, one is left with manifestly 
smooth functions, provided that iV > 2, an obviously acceptable "restriction" 
in our context. This concludes the proof of the smoothness of the coefficients 
B(N,i,v-v v c {v \e ). 

Let us now come to the proof of the statement of Theorem ^ which says 
that the source of a phase transition can only be the second term in square 
brackets in Eq. ffTHj) . which is of topological meaning. To this purpose we have 
to resort to the Main Theorem of paper I and to its Corollary 1. 

Let us assume that only one critical value v c exists in a given interval [vq, Vi]. 
After Sard theorem 171 ]. at any finite N there is a finite number of isolated 
critical points on Sjv« c - For any arbitrarily small 5 > 0, the Main Theorem of 
paper I and its Corollary 1 still apply to the two subintervals [vq,v c — 5} and 
[■y c + 5, Vx]. In order to understand why a breakdown of uniform boundedness in 
iV of \d k S^ /dv k \ f° r k = 3 or k = 4 can be originated only by the topological 
term in r.h.s. of Eq. (fT£|) . we consider each critical point q£ on Ejv« c enclosed 
in a small pseudo-cylindrical neighborhood T(qc \so) of thickness Eq and we 
take e < Eq arbitrarily close to Eq. From Morse theory, we know that passing 
a critical value v c entails that for each critical point of index k{ a /cj-handle 
H N ' ki is attached to M v <Vc so that the following diffeomorphism holds 

M (Vc+£0) « M (Vc . £0) |J H N ' kl |J H N > k > • • • |J H N > k » (49) 

01 02 0n 

where a fcj-handle in N dimensions (0 < ki < N) is the product of two disks, 
a fcj-dimensional disk D ki , and another (N — /cj)-dimensional disk D^ N ~ ki > s.t. 
H N ' hi = D ki x D^ N ~ ki \ and where |J^. stands for the attachment of H N ' ki to 
M(„ c _ £o ) through the embedding fa : § fci_1 x D N ~ ki — > dM^ Vc _ So) (where S is 



an hypersphere; details can be found in 
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3,0 



D. 



The excision of the pseudo-cylindrical neighborhoods T(q l c ,Eo) of all the 
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critical points q£ G ^m c implies that all the manifolds M v := M v \ 
cnt.pts. Y^ l \eo) with A^ — e < v < Nv c + e are free of critical points 
and, consequently, are diffeomorphic. In fact, for any v,v' G M such that 
Nv c — 6q < v < v' < Nv c + Eq, M v is a deformation retraction of M„> through 
the flow associated with the vector field X = -VV N /\\ W N \\ 2 QQ- 

Now, defining M(v,N) = vol(M v ) and T{v,N) = vo\[[j*J rit - pts - ,e )}, 
equation (fTSj) becomes 



v) = ^log[M(v,N) + T(v,N 
1 log \M(v, AO] + ^ log 



1 + 



£(v,N) 
M{v,N) 



(50) 



AT °L v > /J jy 

By applying the Main Theorem of paper I and its Corollary 1 to the first term 
in the r.h.s. of the equation above, we know that j^\d k log [M(u, A 7 ")] /dv k \ for 
k = 1, ... ,4, are uniformly bounded in N, and thus no phase transition can 
be attributed to this term. 

Then, let us consider the second term of the r.h.s. of the equation above. 
By computing its first derivative we obtain 



d 1 

duN 



log 



1 + 



M(v,N) 



i r 



NM+T NM+T 




(51) 



where (M / M) stands for [dM(v ,N)/dv]/M(v,N). After Lemma 1, (M / M) 
is uniformly bounded in A^ and therefore so does the second term in the r.h.s. 
of Eq.ffoTj). Whence, if \dS N ^ /dv\ were to grow with N this could not be due 
to the term ~M(v, N). 

Then we compute the second derivative 

d 2 1 



dv 2 N 



log 



\ r(y,N) 

M(v,N) 



1 T" 1 TM(M /M) + TV 
NM + T + N (M + T) 2 



r 



NM + T 



M \ 1 d M 
W I + Nfa \ W 



(52) 



Again, we can observe that the uniform boundedness with A^ of both (M / M) 
and (d/dv)(M'/M) = (M" fM) - (M' /M) 2 - after Lemma 1 - entails that if 
\d 2 Stf / dv 2 \ were to grow with N, this could not be due to the term M(v, N). 



25 



Similarly, after a lengthy but trivial computation of the third and fourth 
derivatives of the second term in the r.h.s. of Eq.([50|). one finds that 
M(v,N) enters the various terms obtained through the ratio M /M and 
through its derivatives / ' M] with k = 1,2,3, thus, as a consequence 
of Lemma 1, the uniform boundedness in N of [dM(v,N)/dv]/M(v,N) and 
£z{[d~M(v,N)/dv]/M(v,N)} with k = 1,2,3, implies that if \d 3 S ( ~ } /dv 3 \ or 
{frS^/dv^ were to grow with N, this could not be due to the term M(v, N). 

In conclusion, the first term within square brackets in Eq. (jl8j) cannot be at 
the origin of a phase transition, nor can it be the third one, which is the sum 
of smooth functions. Only the second term of the r.h.s. of Eq. ()18|) . which is in 
one-to-one correspondence with topology changes of the M v , can originate an 
unbound growth with N of a derivative \d k / dv k \ for some k, thus entailing 
a phase transition. □ 

Remark 2. A comment about the above considered volume splitting is in order. 
At any finite N , the volume M(v, N) is a smooth function ofv, as is the entropy 
Sj^\v). Since the term ^2f r =Q A(N,i } e )g i jj, i (M v ) entering Eqs. and HJJ[ ) 
depends on the integer valued functions fii(v), one could erroneously think that 
the presence of the functions Hi(v) in this term conflicts with the smoothness 
of volume and entropy. Of course, at finite N, smoothness of volume and 
entropy is not lost. In fact, the term Y2f=o A{N, i, Eo)gi[j J i(M v ) is constant 
in any open interval (vc^ + £q, Vc 1 — £q) and the functions v >— > B(N, i, v) 
smoothly connect in the interval (vc^ —£o, Vc^+Sq) the values that the function 
^2f =0 A(N,i,eo)gi^i{M v ) takes in the intervals (vc 1 + £o, v c^ — £q) and 

(v u c (v) +£ y c {v)+1 -E ). 

Loosely speaking, 'Y^ li A{N : i : £Q)gifii{M v ) + J2 n B(N,i(n),v) is shaped as a 
"staircase" with "rounded corners". 

Remark 3. About the applicability domain of the Main Theorem proved in 
the present paper, note that V is required to be a finite range potential and 
a good Morse function. The former is a physical assumption, the latter a 
mathematical property of V . Finite range potentials are typical in condensed 
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matter systems, where even Coulomb interactions are effective only at a finite 
distance because of the Debye shielding. The mathematical property of being a 
good Morse function is absolutely generic, in fact it requires that a potential 
function is bounded from below and that the Hessian of the potential is non- 
degenerate (i.e. its eigenvalues never vanish). Moreover, given a real-valued 
function f of class C 2 defined on an arbitrary open subset X ofM> N , the map- 
ping x >—> f{x) — {a\X\ + • • • + onXn) '■ X >— > R is nondegenerate for almost all 
(ai, . . . , on) € ^ N (see Chapter 6 of Ref. xlJp] ). This means that nondegeneracy 
is generic whereas degeneracy is exceptional. Continuous symmetries are the 
only physically relevant source of degeneracy, however this kind of degeneracy 
can be removed by adding a generic term {a\X\ + • • • + o^xn) to the potential 
with an arbitrarily small vector (ai, . . . , a/v). This removal of degeneracy is 
a rephrasing, within the framework of Morse theory, of a standard procedure 
undertaken in statistical mechanics to explicitly break a continuous symmetry, 
that is the addition of an external field whose limit to zero is taken after the 
limit N — > oo. 

Remark 4. Let us briefly compare meaning and strength of the Main Theorems 
of papers I and II, which we denote by MT-I and MT-II respectively. The 
proof of MT-I is preliminary to, and independent of, the proof of MT-II. On 
the other hand, it is important to note that the hypotheses under which MT-II 
applies are fulfilled by a very broad class of physically meaningful potentials, 
at variance with the case of MT-I. In fact, MT-I applies in presence of a 
potential energy density interval of finite length which is definitevely (that is 
asymptotically in N) free of critical points, a rather restrictive assumption 
which we could hardly relax in the present demonstration scheme. On the 
contrary, we can approach the problem from a complementary point of view 
allowing the existence of critical points of the potential - as is assumed in the 
hypotheses of MT-II - with the only limitation to a linear growth with N of 
the number of critical values of V in a potential energy density interval of 
finite length. This assumption is suggested by what happens in lattice systems 
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where critical levels are separated by a finite minimum energy amount, as is 
the case of "spin flips" in the ld-XY model mean- field XY model 



p-trig model L r ll], or of elementary configurational changes that, in lattices and 
fluids, correspond to the appearance of a new critical value ofV at an energetic 
cost which is independent of N. 

It would be impossible to prove MT-II without resorting to MT-I, though 
MT-II is much stronger than MT-I. Actually, the link established by MT-II in 
Eq. hltfy . between thermodynamic entropy and topology, tells that only through 



suitable variations with N of the v-patterns of^2 i=0 A(N,i,£o)gi/ii(M- l 
of the derivatives \d k S^ N /dv k \ can cease to be uniformly bounded in N from 
above, thus giving rise to a phase transition. This is an important hint for a 
future rigorous investigation of the sufficiency conditions for both MT-I and 
MT-II. The combination of MT-I and MT-II provides clear evidence of the 
relevance of topology for the phenomenon of phase transitions. 

4. FINAL REMARKS 

Let us conclude with a few general comments. Earlier attempts at introduc- 
ing topological concepts in statistical mechanics concentrated on macroscopic 
low-dimensional parameter spaces. Actually this happened after Thorn's re- 
mark that the critical point shown by the van der Waals equation corresponds 



some 



to the Riemann-Hugoniot catastrophe 



22j. Hence some applications of the 



theory of singularities of differentiable maps to the study of phase transitions 



followed 



23|. An elegant formulation of phase transitions as due to a topologi- 



cal change of some abstract manifold of macroscopic variables was obtained by 



251 ] and deserves special attention 



using the Atiyah-Singer index theorem 
because it applies to the 2d Ising model, whose phase transition is associated 
with a jump of the Atiyah index of some suitable vector bundle. This shows 
that also for discrete variables systems, like spin systems, topological concepts 
can be useful in the study of phase transitions, provided that the relevant 
manifolds are identified. 
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The Main Theorem, that we have proved above, makes a new kind of link 
between the study of phase transitions and differential topology. In fact, in 
the present work we deal with the high-dimensional microscopic configuration 
space of a physical system. The level sets of the microscopic interaction po- 
tential among the particles - or the regions of configuration space bounded 
by them - are the configuration space submanifolds that necessarily have to 
change their topology in correspondence with a phase transition point. The 
topology changes implied here are those described within the framework of 
Morse theory through attachment of handles jl8|. 

We have explicitly investigated these topology changes in some particular 
models. The results so far obtained, already reported in the literature, are: 
i) numerical results on the lattice if 4 model |?|, ii) exact analytical results 
on the two models considered in Refs. hJ Q], Hi) analytical results on the 
coupled rotators model. Although the models in item (ii) do not fulfil the 
condition of short-range interactions assumed by MT-I and MT-II in their 
present formulations, we already get a coherent scenario illustrating in practice 
how topology changes in presence and in absence of phase transitions. In 
fact, we observe that when the ^-pattern of the Euler characteristic x{v) ~ 
which is the probe that we use to detect the variations of topology with v - 
approaches a smooth curve as N increases, then phase transitions are absent; 
this is the case of x(S^) for the one- dimensional (p A model and of x(M v ) for 
the one-dimensional chain of coupled rotators. At variance, sharp jumps or 
"cuspy" w-patterns of x(£v) or °f x{M v ) are associated with phase transitions, 
in qualitative agreement with what is expected after MT-II; in Ref . [h | it is also 
shown that first and second order phase transitions are signaled by markedly 
different w-patterns of x{M v ). 

Notice that in our approach the role of the potential V is twofold: it deter- 
mines the relevant submanifolds of configuration space and it is a good Morse 
function on the same space. However, for example, in the case of entropy 
driven phase transitions occurring in hard sphere gases, the fact that the (sin- 
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gular) interaction potential cannot play any longer the role of Morse function 
does not mean that the connection between topology and phase transitions is 
lost, it rather means that other Morse functions are to be used. Just to give 
an idea of what a good Morse function could be in this case, let us think of 
the sum of all the pairwise euclidean distances between the hard spheres of 
a system: it is real valued, it has a minimum when the density is maximum, 
that is for close packing, meaning that this function is bounded below. The 
discussion of non- degeneracy is more involved and here would be out of place, 
let us simply remark that Morse functions are dense and degeneracy is easily 
removed when necessary. 

The topology of configuration space submanifolds makes also a subtle link 
between dynamics and thermodynamics because it affects both of them, the 
former because it can be seen as thegeodesic flow of a suitable Riemannian 
metric endowing configuration space [9j, the latter because an analytic (though 
approximate) relation between thermodynamic entropy and Morse indexes of 
the critical points of configuration space submanifolds can be worked out . 

Moreover, there are "exotic" kinds of transitional phenomena in statisti- 
cal physics, like the glassy transition of amorphous systems to a supercooled 
liquid regime, or the folding transitions in polymers and proteins, which are 
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27| which 



qualitatively unified through the so-called landscape paradigm 
is based on the idea that the relevant physics of these systems can be under- 
stood through the study of the properties of the potential energy hypersurfaces 
and, in particular, of their stationary points, usually called "saddles". That 
this landscape paradigm naturally goes toward a link with Morse theory and 
topology has been hitherto overlooked. However, though at present our Main 
Theorem only applies to first and second order phase transitions, the topo- 
logical approach seems to have the potentiality of unifying the mathematical 
description of very different kinds of phase transitions. 



30 



5. ACKNOWLEDGMENTS 



The authors wish to thank L. Casetti and S. Schreiber for comments and 
suggestions. A particularly warm acknowledgment is addressed to G. Vezzosi 
for his continuous interest in our work and for many helpful discussions and 
suggestions. 



[1] L. Onsager, Phys. Rev. 65, (1944) 117. 

[2] C.N. Yang and T.D. Lee, Phys. Rev.87, (1952) 404. 

[3] D. Ruelle, Thermodynamic formalism, Encyclopaedia of Mathematics and its Applications, 

(Addison- Wesley, New York, 1978). 
[4] H.O. Georgii, Gibbs Measures and Phase Transitions, (Walter de Gruyter, Berlin, 1988). 
[5] L. Caiani, L. Casetti, C. Clementi and M. Pettini, Phys. Rev. Lett. 79, (1997) 4361. 
[6] R. Franzosi, L. Casetti, L.Spinelli and M. Pettini, Phys. Rev. E60, (1999) R5009. 
[7] R. Franzosi, M. Pettini, and L.Spinelli, Phys. Rev. Lett. 84, (2000) 2774. 
[8] L. Casetti, E.GD. Cohen and M. Pettini, Phys. Rev. Lett. 82, (1999) 4160. 
[9] L. Casetti, M. Pettini, and E.GD. Cohen, Phys. Rep. 337, (2000) 237-341. 
[10] L. Casetti, M. Pettini, and E.GD. Cohen, J. Stat. Phys. Ill, (2003) 1091. 
[11] L. Angelani, L. Casetti, M. Pettini, G. Ruocco, and F. Zamponi, Europhys. Lett. 62, (2003) 

775; Phys. Rev. E71, (2005) 036152. 
[12] R. Franzosi, M. Pettini, and L.Spinelli, Topology and Phase Transitions I. Preliminary Results, 

archived in: |math-ph /0505057. 
[13] R. Franzosi, and M. Pettini, Phys. Rev. Lett. 92, (2004) 060601. 

[14] M. Morse and S. S. Cairns, Critical Point Theory in Global Analysis and Differential Topology, 

(Academic Press, New York, 1969). 
[15] H. Federer, Geometric Measure Theory, (Springer, New York 1969), p. 249. 
[16] D. Ruelle, Statistical Mechanics. Rigorous results, (Benjamin, Reading, 1969). 
[17] R.S. Palais and C. Terng, Critical Point Theory and Submanifold Geometry, (Springer, New 

York 1988). 

[18] MW. Hirsch, Differential Topology, (Springer, New York 1976). 



31 



[19] J. Milnor, Morse Theory, (Princeton University Press, Princeton, 1973). 
[20] P. Laurence, ZAMP 40, (1989) 258. 

[21] R. Bott and J. Mather, Topics in Topology and Differential Geometry, in Battelle Rencontres, 
Eds. CM. De Witt and J. A. Wheeler, p. 460. of the Main Theorem of paper I. The combination 
of the Main Theorem of paper I and Theorem ^ of the present paper 

[22] R. Thorn, in Statistical Mechanics, Eds. S.A. Rice, K.F. Freed, and J.C. Light, (University of 
Chicago Press, 1972), p.93. 

[23] T. Poston and I. Stewart, Catastrophe Theory and its Applications, (Pitman Press, London, 
1978), and references therein quoted. 

[24] M. Rasetti, Topological concepts in the theory of phase transitions, in Differential Geometric 
Methods in Mathematical Physics, Ed. H.D. Dobner (Springer- Verlag, New York, 1979). 

[25] M. Rasetti, Structural Stability in Statistical Mechanics, in Springer Tracts in Math., Ed. W. 
Giittinger (Springer- Verlag, New York, 1979). 

[26] F.H. Stillinger, Science 267, (1995) 1935. 

[27] S. Sastry, P.G. Debenedetti and F.H. Stillinger, Nature 393, (1998) 554. 



32 



